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Abstract 

In this approach to discrete quantum gravity the basic structural 
element is a covariant causal set (c-causet). The geometry of a c- 
causet is described by a shell-sequence that determines the discrete 
gravity of a universe. In this growth model, universes evolve in dis¬ 
crete time by adding new vertices to their generating c-causet. We first 
describe an inflationary period that is common to all universes. After 
this very brief cycle, the model enters a multiverse period in which the 
system diverges in various ways forming paths of c-causets. At the be¬ 
ginning of the multiverse period, the structure of a four-dimensional 
discrete manifold emerges and quantum mechanics enters the picture. 
A natural Hilbert space is defined and a discrete, free Dirac opera¬ 
tor is introduced. We determine the eigenvalues and eigenvectors of 
this operator. Finally, we propose values for coupling constants that 
determine multiverse probabilities. These probabilities predict the 
dominance of pulsating universes. 


1 Causal Sets 

This section presents a brief review of causal sets (causets). For more back¬ 
ground and details, we refer the reader to mmn- We call a finite partially 
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ordered set {x,<) a causet and interpret a < 6 in x to mean that b is in the 
causal future of a ilElEI. A labeling for a causet x of cardinality |x| is a 
bijection 

X -)■ { 1 , 2 ,..., |x|} 

such that a,b E X with a < b implies i{a) < i{b). A labeling for x can be 
thought of as a “birth order” for the vertices of x. Two labeled causets x, y 
are isomorphic if there is a bijection (j): x ^ y such that a < 6 in x if and 
only if 0 (a) < 0 ( 6 ) in y and i [ 0 (a)] = i{a) for all a G x. A causet is covariant 
if it has a unique labeling (up to isomorphism) and we call a covariant causet 
a c-causet M- Denote the set of c-causets with cardinality n by V'^ and the 
set of all c-causets by V. It is shown in HH that the cardinality |P^| = 2"' 
n = 1 , 2 ,.... Two vertices a,b E x are comparable if a < 6 or 6 < a. The 
height h{a) of a E x is the cardinality, minus one, of a longest path in x that 
ends with a. It is shown in | 2 ] that a causet x is covariant if and only if 
a,b E X are comparable whenever h{a) ^ 6.(6). 

If X G P' we call the sets 

Sj{x) = (a G x: 6(a) = j} , j = 0 , 1 , 2 ,... 

shells and the sequence of integers Sj{x) = IS'j^x)!, j = 0,1,..., is the shell 
sequence for x. A c-causet is uniquely determined by its shell sequence and 
we think of {si(x)} as describing the “shape” or geometry of x. In this model, 
we view a c-causet as a framework or scaffolding for a possible universe at 
a hxed time. The vertices of x represent small cells that can be empty or 
occupied by a particle. We shall later consider a growth model in which 
universes evolve in time by adding new vertices. 

Let X = (ai, a2,..., an} E V'n, where the subscript i of ai is the vertex 
label, f'(aj) = i. A path is a sequence 


y ajj^ajj ■ ■ ■ 


in X starting at a^^ and moving along successive shells until a^^ is reached. 
We dehne the length of 7 by [T] 




D-i) 

-i=2 


For a,b E X with a < b, a geodesic from a to 6 is a path from a to b that has 
the shortest length [T]. Clearly, if a < 6, then there is at least one geodesic 
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from a to b. If a, 6 G x are comparable and a < b say, then the distance 
d{a,b) is the length of a geodesic from a to 6 pQ. It is shown in PQ that if 
a < c < b, then d{a, b) < d{a, c) + d{c, b). Thus, the triangle inequality holds 
when applicable so d{a, b) has the most important property of a metric. It is 
also shown in [1] that a subpath of a geodesic is a geodesic. 

For conciseness, let us refer to a vertex by its label. If there are j geodesics 
from vertex 1 to vertex n, we dehne the curvature K{n) at n to be K{n) = 
i — 1 [I]. One might argue that the curvature should be a local property 
and should not depend so heavily on vertex 1 which could be a considerable 
distance away. However, if there are a lot of geodesics from 1 to n, then by 
the last sentence of the previous paragraph, there are also many geodesics 
from other vertices to n. Thus, the dehnition of curvature is not as dependent 
on the initial vertex 1 as it hrst appears. Assuming that particles tend to 
move along geodesics, we see that K[n] gives a measure of the tendency for 
vertex n to be occupied. In this way, the geometry oi x eV determines the 
gravity at each vertex of x. 


2 Inflationary Period 


In this section and the next we shall refer to a vertex in a c-causet by its label. 
Thus, for X E V'^ we have x = {1,2,... ,n}. According to recent observa¬ 
tions and cosmological studies, immediately after the big bang the universe 
expanded exponentially during an inflationary period. We can describe this 
period by the c-causet Xq^q with shell sequence (1, 2,4, 8 ,..., 2"'). We see that 
Xo,o has X]r=o ~ vertices. We can obtain an approximate value for n 
by noting that the number of photons in the universe is approximately 10®®. 
We assume that most of these photons were created during the inflationary 
period and that they still exist as the cosmic microwave background. We 
also assume that during the inflationary period, the only particles that ex¬ 
isted were photons and possibly “dark photons” which result in dark energy. 
Moreover, including empty cells through which photons can move we esti¬ 
mate that the number of vertices produced during the inflationary period to 
be about 10®®. Hence, 2"+® ^ 10®® so that 


InlO 

n = 93 --1 

ln2 


308 


Postulating that each shell is hlled in one Planck instant given by about 10 ^® 
seconds, we conclude that the inflationary period lasted only about 3 x lO”"^® 
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seconds. 

It has also been experimentally verified that the universe, in the large, 
is quite fiat and homogeneous. This suggests that during the infiationary 
period, the curvature should be essentially zero. We now indicate why this is 
true for our present model by considering some small cases. As an example of 
a toy inflationary universe, suppose x has shell sequence (1, 2,4, 8,16). With 
shells delineated by semi-colons, we label the vertices as follows: 

(1; 2, 3; 4, 5, 6, 7; 8, 9,10,11,12,13,14,15; 16,17,..., 31) 

Then (1,2) and (1,3) are geodesics, d{l,2) = 1, <7(1,3) = 2 and K(2) = 
K{3) = 0. Similarly, (1,2,4), (1,3,4), (1,3,5), (1,3,6) and (1,3,7) are 
geodesics, <7(1,4) = \/5, <7(1,5) = \/8, <7(1,6) = vTS, <7(1,7) = \/^ and 
77(4) = 1, k{5) = 77(6) = 77(7) = 0. Table 1 displays the distances squared 
and curvatures for the vertices of x. 


Vertex 7 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

<7(1,7)^ 

1 

4 

5 

8 

13 

20 

17 

22 

29 

39 

45 

56 

69 

84 

61 

70 

81 

K(i) 



1 




1 


1 






1 


1 


Vertex 7 

19 


21 

22 

23 

24 

25 

26 

27 

28 

29 


31 

<7(1,7)^ 

92 

105 

118 

133 

148 

165 

184 

205 

228 

253 

280 

309 

340 

K(i) 

0 

1 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 


Table 1 (Distances Squared and Curvatures) 


Continuing this example to the next shell with 32 vertices, we find that 
77(7) =0 except for i = 32,34,36,38,40,42,44,46 in which case 77(7) = 1. 
We conclude that in this case, the curvature is zero except for a sprinkling 
of ones which are about one-fourth the total. 

That this conclusion always holds is reinforced by considering the geodesics 
for the c-causet with shell sequence (1,2,4,8,16,32). We denote a geodesic 
with vertices having labels 7i, 72 ,..., 7^ by (7i, 72 ,..., in)- The geodesics ter¬ 
minating at vertices in shell 1 are: (1,2), (1,3). Those terminating at vertices 
in shell 2 are: (1,2,4), (1,3,4), (1,3,5), (1,3,6), (1,3,7). 
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For shell 3 we have: (1,3,5,8), (1,3,6,8), (1,3,6,9), (1,3,6,10), (1,3,7,10), (1,3,7,11), 
(1,3,7,12), (1,3,7,13), (1,3,7,14), (1,3,7,15). 

For shell 4 we have: (1,3,7,12,16), (1,3,7,13,16), (1,3,7,12,17), (1,3,7,12,18), 
(1,3,7,13,18), (1,3,5,13,19), (1,3,7,13,20), (1,3,7,14,20), (1,3,7,14,21), (1,3,7,14,22), 
(1,3,7,15,22), (1,3,7,15,23), (1,3,7,15,24),.(1,3,7,15,31). 

For shell 5 we have: (1-15,23,32), (1-15,24,32), (1-15,24,33), (1-15,24,34), 
(1-15,25,34), (1-15,25,35), (1-15,25,36), (1-15,26,36), (1-15,26,37), (1-15,26,38), 
(1-15,27,38), (1-15,27,39), (1-15,27,40), (1-15,28,40), (1-15,28,41), (1-15,28,42), 
(1-15,29,42), (1-15,29,43), (1-15,29,44), (1-15,30,44), (1-15,31,48), (1-15,31,49), 
..., (1-15,31,63). 


This pattern continues indehnitely and we conclude that the curvatures are 
all zero or one with approximately 3/4 zeros and 1/4 ones. 

Let m, m+1 ,..., m+n be the labels of the vertices in a shell for a c-causet. 
We call m + n the endpoint of the shell and \m + the midpoint where [rJ 
denotes the integer part of r. Let x be an exponentially growing c-causet 
with shell sequence (1, 2,,,, 4,..., 2"'). It follows by induction on the shell 
number that every geodesic (ii,i 2 , • • •*„), n > 3, of a; has ii,i 2 , • • • ,*n -2 as 
endpoints. Also, in-i is an endpoint if in is at least as large as the midpoint 
of its shell and moreover, in-i is at least as large as the midpoint of its shell. 
This gives a fairly complete description of the geodesics in x. 

We have seen that exponentially growing c-causets have a uniformly low 
curvature. It is of interest that this is a characteristic property of exponen¬ 
tial growth. That is, a c-causet that grows at a slower pace does not have 
this property. To illustrate this, let y be the c-causet with shell sequence 
(1, 2, 3,..., 10) so y is growing very slowly. Notice that y has 55 vertices so y 
is about the same size as the exponentially growing c-causet considered pre¬ 
viously. However, the curvatures are considerably different as the following 
table shows. 


Vertex i 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

m 

0 

0 

1 

0 

0 

0 

1 

0 

0 

2 

0 

1 

0 

0 

2 

2 


Vertex i 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

m 

0 

1 

0 

0 

3 

2 

1 

0 

0 

0 

0 

2 

0 

1 


5 



Vertex i 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

m 

1 

0 

1 

0 

0 

3 

2 

0 

1 

2 

0 

0 

0 


Vertex i 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

m 

0 

2 

4 

3 

0 

0 

1 

0 

0 

0 

0 


Table 2 (Curvatures in c-causet y) 


The inflationary period lasts until the system arrives at the c-causet a;o,o 
with shell sequence (1, 2,4,..., 2"') with n ~ 308. For this n the system can 
no longer sustain the immense energies, pressures and densities. The infla¬ 
tionary period terminates and the system goes through a phase transition. 
At this point most of the vertices (cells) of a universe have been produced 
and the exponential growth ceases. The system resembles a huge sea of ver¬ 
tices and the vertex production continues above this sea at a much slower 
pace. The model then enters what we call the multiverse period. This is the 
period in which we humans now And ourselves. 

3 Multiverse Period 

The multiverse period begins with the c-causet a;o,o- Then the system diverges 
in various ways forming different c-causet paths starting at the common c- 
causet a;o,o- Each path corresponds to a distinct universe history and these 
evolving universes form a multiverse model. During the inflationary period, 
there are no physical principles operating except the geometry given by the 
shell sequence. In particular, the curvature is essentially zero and there is 
no concept of dimension. As we shall see, at the beginning of the multi¬ 
verse period, the structure of a 4-dimensional discrete manifold emerges and 
quantum mechanics enters the picture. In this sense, the much sought after 
“theory of everything” is quantum mechanics itself. Moreover, very early in 
the multiverse period, high curvatures are created which explains why matter 
collects in certain places and not in others. However, these curvatures are 
local and since most of the vertices of any of the universes have already been 
formed during the inflationary period, the universes remain relatively flat. 
This is because once a curvature has been established it remains unchanged 
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at later times. In the rest of this section and the next, we shall make these 
statements precise. 

An element a of a causet x is maximal if there is no 6 G a; with a < b. 
Thus, a is maximal in x if there are no elements of x that are in the causal 
future of a. li x,y G V' we say that x produces y ii y \s obtained by hrst 
adjoining a maximal element a to a; and then adjoining a second maximal 
element 6 to a;U {a} so that y = xU {a, h} where a ^ x, h ^ x. If x produces 
y, we write x ^ y and say that y is an offspring of x It is easy to see that any 
X G P' has precisely four distinct offspring. In fact, if x has shell sequence 
(so, si,..., Sn), then the offspring of x have shell sequences: 

(sq. Si, ..., s„ + 2) 

(sq. Si, ..., s„ + 1 ,1) 

(■Sqj Si) ■ ■ ■) Sn, 2) (3-1) 

(SQ) Si, . . . , Sn, 1, 1) 

We number the c-causets in the multiverse period recursively as follows. The 
four offspring of Xo,o in the order given by fl3.ip are denoted by Xi^O) 2^i,i) 2 ^ 1 , 2 ) 
Given c-causet Xnj the offspring in the order given by (13.1 p are Xn+i, 4 j+fc, 
k = 0,1, 2, 3. For example, the offspring of Xi^O) Xi^i,Xi^ 2 i 2 ^ 1,3 are: 


2^2,0) 2^2,1) 2^2,2) 2:2,3 
2^2,4) 2^2,5) 2:2,6) 2:2,7 
2 ^ 2 , 8 ) 2 : 2 , 9 ) 2 : 2 , 10 ) 2 : 2,11 

2:2,12) 2 : 2 , 13 ) 2 : 2 , 14 ) 2:2,15 


respectively. 

We use the notation Vo = { 2 : 0 , 0 } ) Pi = { 2 : 1 , 0 ) 2 : 1 , 1 ) xi, 2 ) xi, 3 } ,..., 

Vn = {xn,j:j = 0,l,...,A^-l} 

We thus see that |Pn| = d"'. We denote the set of all c-causets in the mul¬ 
tiverse period by P = UfLQVn- We can view V as & 4-dimensional discrete 
manifold in which the “tangent vectors” at Xn,j consist of the four edges 

^n,j 2:n+l,4j-|-A:)) k 0, 1, 2, 3 

We think of dnj, k = 0,1, 2, 3, as the four directions originating at x„j. 
For simplicity, we consider xo,o as the c-causet with a single vertex that we 
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denote by 0. Thus, we identify a;o,o with the c-causet {0}. In this way w 
have |a;o,o| = 1 and the shell sequence of Xq^ is simply (1). In a similar 
way |a;ij| = 3, j = 0,1,2,3 and the shell sequence of xij, j = 0,1,2,3 are 
(3), (2,1), (1, 2), (1,1,1), respectively. In general \xn,j\ = 2n+l and the shell 
sequences of Xnj, = 0,1,...,4"' — 1 are 

(2n + 1), (2n, 1), (2n - 1, 2), (2n - 1,1,1),..., (1, !,...!) 
respectively. 

For the multiverse period, we dehne paths, distances, geodesics and cur¬ 
vature as before except now we begin with the vertex 0. To show how large 
curvatures can be generated, we consider an example of an ideal pulsating 
universe Xp. Although our universe may not be described exactly by Xp, we 
are proposing that it may be similar to Xp. The c-causet Xp has shell sequence 

(1,2,4, 2,4, 6,4, 2,4, 6, 8, 6,4,2,..., 6,4, 2) (3.2) 

As before, we denote the vertices by their labels. The vertices are listed as 
follows where we stop at vertex 72: 


(0; 1, 2; 3,4, 5, 6; 7, 8; 9,10,11,12; ■ ■ ■ ; 65, 66, 67, 68, 69, 70, 71, 72) 


Examples of two geodesics are (0,2,4, 7) and (0,2,5, 7). It follows that 
(i(0, 7) = vTf and K{7) = 1. Also, (0,2, 5, 8) is a geodesic, d{0,8) = \pn 
and k{8) = 0. Table 3 displays the distances squared and curvatures for the 
vertices of Xp up to 72. 


Vertex i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

d{0, i)^ 

1 

4 

5 

8 

13 

20 

17 

22 

21 

26 

31 

38 

35 

40 

47 




1 




1 


1 

2 



2 


1 


Vertex i 

16 

17 

18 

19 


21 

22 

23 

24 

25 

26 


28 

29 

(i(0, i)^ 

54 

63 

74 

63 

70 

79 

88 

79 

86 

83 

88 

95 

102 

97 

K(t) 

0 

0 

0 

2 

0 

1 

0 

3 

0 

3 

3 

4 

0 

3 
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Vertex i 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

d{0, 

104 

111 

118 

127 

138 

127 

134 

143 

152 

163 

174 

K(,t) 

8 

5 

0 

0 

0 

6 

0 

1 

0 

1 

0 

Vertex i 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

d{0, 

187 

202 

177 

188 

199 

210 

223 

236 

213 

224 

235 

m 

0 

0 

0 

1 

1 

0 

1 

0 

2 

3 

2 

Vertex i 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

(i(0, iY 

246 

229 

238 

233 

238 

245 

254 

247 

254 

261 

270 

m 

0 

2 

2 

2 

2 

2 

5 

2 

5 

2 

8 

Vertex i 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 


d{0, iY 

279 

290 

277 

286 

295 

304 

315 

326 

339 

354 


K(i) 

5 

11 

2 

11 

14 

5 

17 

11 

11 

11 



Table 3 (Distances Squared and Curvatures for Xp) 


We call the set of vertices in Xp between 2s in the shell sequence fl3.2p a 
pulse. Starting with the first 2 in (I3.2p . let po,pi,p 2 ,... be the pulses. That 
is, po = {1,2,...,7,8}, Pi = {7,8,..., 23, 24}, p 2 = (23, 24,..., 53, 54} ,.... 
We see that |po| = 8, |pi| = 18, |p 2 | = 32, |p 3 | = 50,.... In general we have 
that 

\Pn\ = 8 + lOn + 4 [1 + 2 +-h (n — 1)] = 8 + lOn + 2(n — l)n = 2{n + 2)^ 

As previously discussed, during the multiverse period, we are assuming that 
two vertices are created every Planck instant. It follows that the time dura¬ 
tion of the nth pulse in Planck instances is 

tn = I \Pn\ - 1 = (n -h 2)^ -h 1 = (n + l)(n 3) 

Thus, to = 3,ti = 8, ^2 = 15, ^3 = 24,.... The time in Planck instants at 
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which the nth pulse begins is computed as follows: 
hj = 1) = 1 + ^0 = 4, T 2 = 1 + to + ti = 12, T 3 = 1 + to + ti + t 2 = 27,... 


In general, we have that 


n—1 


n—1 


n—1 


Tn = 


1 + 'y ^ ti — 1 + y + i)(t + 3 ) — 1 + y ^ (3 + 4 t + 


i=0 


i=0 


i=0 


= 1 + 3n + 2(n — l)n + 


(n — l)n( 2 n — 1 ) 
6 


= 1 + n( 2 n + 1 ) + 


(n — l)n( 2 n — 1 ) 
6 


The time of our present is approximately 10®° Planck instants. If we 
are in the nth pulse now, then Tn ~ 10®°. Since n is large, we have that 
Tn ~ n°/3. Hence, n°/3 ~ 10®° so that n ~ 3^t° x 10^°. We conclude 
that the universe Xp is about in its 10^°th pulse. The duration of this pulse is 
tn ~ 10*^° Planck instants which is approximately 10“° seconds. If this is near 
the behavior of our own universe, it would be close to being imperceptible. 
Finally, the number of new vertices produced in each universe during the 
multiverse period is about 2 x 10®° This is much smaller than the number 
of vertices produced during the inflationary period. It has been estimated 
that the number of massive particles in our universe is about 10®°. Some 
of these particles were produced during the multiverse period, but most of 
them came from photons produced during the inflationary period. 


4 Quantum Mechanics Emerges 

This section discusses the emergence of quantum mechanics at the beginning 
of the multiverse period. The set V of odd cardinality c-causets together with 
the production relation gives a tree (P, ^•). We have seen in Section 3 that 
we can denote the c-causets in V by Xn,j-, n = 0,1, 2,..., j = 0,1,..., 4*^ — 1. 
Moreover, we have four directions (tangent vectors) 

^n,j ^n+l,4:j+k) j k 0, 1, 2, 3 

emenating from each node (c-causet) Xn,j- We say that two tangent vectors 
are incident if they have the form (x, y) and (y, z). We also use the notation 
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Each X E V except a:o,o has a unique producer and each x E V has exactly 
four offspring. In particular Xn,j Xn+i,Aj+ki k = 0,1, 2, 3 and we interpret 
the tree (P, -^) as a sequential growth process. 

An n-path in P is a sequence oj = uqUi ■ ■ - Un where Ui E Vi and Ui -E 
cji+i. We denote the set of n-paths by fin and interpret an a; G fi as the 
history of the universe Un at which u terminates. We can also consider an 
n-path as a sequence of tangent vectors 


u = 


“0,0 “iji 


■■■d 


^n — 1 

n-lJri-1 


where each tangent vector is incident to the next. Since each Xnj E Vn has 
a unique history, we can identify Vn with fin and we write Vn ^ kin- 

We now describe the evolution of the multiverse in terms of a quantum 
sequential growth process. In such a process, the probabilities of compet¬ 
ing geometries are determined by quantum amplitudes. These amplitudes 
provide interferences that are characteristic of quantum systems. A tran¬ 
sition amplitude is a map a: V x V ^ C satisfying a{x,y) = 0 if x y 
and^^g.p a(x, y) = I for all x eV. We call a a unitary transition amplitude 
(uta) if a also satishes = 1 for all x eV. We conclude that 

a uta satishes 


3 

^ ^ ^n+l,4j+fc) ^ ^ y ^n+l,4j+k)\ 1 


A:=0 


fc =0 


(4.1) 


for all u = 0,1, 2,..., j = 0,1,..., d"" — 1. Using the obvious notation, we 
can also write fld.ip as 


E“(A)=Ei2«i)r=i 

A:=0 A:=0 

It is shown in |3lll] how utas can be constructed. We call 

^n,j ^n+l,4:j+k) y k 0) 2, 3 

the coupling constants for a. As of now, the coupling constants are undeter¬ 
mined. It is hoped that their values can be found so that general relativity 
becomes an approximation to the present theory. We have thus obtained a 
kind of quantum Markov chain in which transition probabilities are replaced 
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by utas. If a is a uta and oj = ojqUJi ■ ■ - oJn ^ we define the amplitude of 
uj to be 

a{oj) = a{ojQ,ui)a{ui,U2) ■ ■ ■a{oJn-i,uJn) 

Also, we define the amplitude of x E Vn to be a{oj) where uj is the nniqne 
path in that terminates at x. 

Let Hn be the Hilbert space 

= L^iQn) = L^iVn) = {f:Vn^C} 

with the standard inner prodnct 

(/,£/) = /(^)^(^) 
xe.Vn 

Let Xn,j be the nnit vector in given by the characteristic fnnction Xx^ ■■ 
Thns, dim. Hn = 4" and {xn,j : J = 0,1,..., 4” — 1} becomes an orthonormal 
basis for Hn- For the remainder of this section, a is a nta with corresponding 
conpling constants We now describe the qnantnm dynamics generated 
by a. Dehne the operators Un '■ Hn —t Hn+i by 

3 

Un^nJ ^ ^ Cn jXn-\-lXj+k 
k=0 

and extend Un to Hn by linearity. It is shown in [3lll] that Un is a partial 
isometry with U^Un = U where the adjoint U* ■ Hn+i —)■ Hn of Un is given 
by 

Un^n+l^ij+k ^n,j^n,ji ^ 0, 1, 2, 3 

As nsnal, a state on Hn is a positive operator p on Hn with tr(p) = 1. A 
stochastic state on Hn is a state p that satishes (pin, In) = 1 where = XvU 
that is, In(a^) = 1 for all x EVn- It is shown in |3] that if p is a state on Hn, 
then UnpU* is a state on Hn+i- Moreover, if p is a stochastic state on Hn, 
then UnpU* is a stochastic state on Hn+i- The qnantnm dynamics is then 
given by the map p i—)■ UnpU*- For farther details, we refer the reader to [1]. 

Althongh Un ■ Hn —t Ffn+i is a partial isometry, it cannot be nnitary be- 
canse Hn and Ffn+i have different dimensions. However, when an additional 
condition is imposed, we can constrnct a related nnitary operator. We say 
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that a and are strong if the matrices given by 


K = 


c° • 

ri,j 

■ 

n,j 

■ 

n,] 

. 

n,j 

■ 

n,j 

C° • 
n,j 

■ 

n,j 

• 

ri,j 



c° . 

'^n,j 

■ 

■ 

■ 

C° • 


(4.2) 


are nnitary. Notice that if a is a nta, then Y^k=o knjl = 1 so Kn,j is nnitary 
if and only if 

+ cljclj) = Re(c° + cljclj) 

= Re(c°_jc3. + = 0 (4.3) 

Thns, fl4.3p is a necessary and snfficient condition for a nta to be strong. 
A nta already satishes a condition similar to fl4.3p bnt not qnite as strong. 
Since J2k=o ^n,j = 1 we have 


1 = 




k=0 

It follows that 






k=0 k'=0 


k=0 


k,k'^0 

k<k' 


Y1 = 0 


(4.4) 


/c,fc'=0 

k<k' 


Notice that fl4.3l) implies fl4.4p . It follows that if Kn,j is nnitary, then a is 
essentially a nta. In fact when 


EKil 


fc =0 


fc =0 




= 1 


If we mnltiply the by a hxed phase factor then becomes 

a nta. 

A direct verihcation shows that the eigenvalnes of Knj are 

lO _ ^0 I I 2 I 3 


-X- A- -A _ 1 

' ^n,j ' ' ^n,j 

\1 ^^0 _ 1 I 2 _ 3 

^n+l,j ^n,j ^n,j ^n,j ^n,j 
\2 = ^0 , 1 _ 2 _ 3 

\3 = pO _ 1 _ 2 , 3 


(4.5) 
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with corresponding unit eigenvectors 


i(l, 1,1,1), |(1, -1,1, -1), i(l, 1, -1, -1), |(1, -1, -1,1), 
respectively. 

The Hilbert space ifn+i can be decomposed into the direct sum 


Hn+l — Q © i^n+1,1 © ' ' ' © i^n+l,4"-l (4-6) 

where dimif„+ij = 4 and an orthonormal basis for n = 0,1, 2 ..is 

{^n+l,4j+k ■ k 0, 1, 2, 3} 

We dehne a stochastic unitary operator W+ij : Hn+ij —)■ with matrix 

representation given by fl4.2p . To be explicit, we have 

hn+l,j*^n+l,4j ^n,j^n-+l,4j+2 T Wj'^n+l,4j+3 

hn+lj^n+l,4j+l W,j^»^+l,4j T WjXn+l,4j+l T Wj3^n+l,4j+2 T Wj'^n+l,4j+3 

hn+lj3^n+l,4j+2 ^n,j^n+l,4:j ^n,j^n+l,4:j+l T C^j^n+l,4j+2 T Wj'^n+l,4j+3 

hn+ljirn+l,4j+3 Wj^«-+l,4i T Wj'2^n+l,4i+l T Wj2^n+l,4i+2 T j’2^n+l,4j+3 


We conclude from our previous work that the eigenvalues of Vn+ij are ^n+ij^ 
fc = 0,1, 2, 3 given by fl4.5p with corresponding unit eigenvectors 


a; 


n+lj' 2(^n+l,4i T lCn+l,4i+l T ^n+l,4i+2 T Xn+l,4i+3) 

a^n+lj' 2 (^’^+li© ^n+l,4j+l T a^n+l,4j+2 ^n+l,4j+3) 
a^n+lj "h ^n+l,4i+l a^n+l,4i+2 a^n+l,4j+3) 


x; 


^n+lj 2(^”+1 i4j ^n+l,4j+l ^n+l,4j+2 T Xn+lj4j+3) 

Finally, we define the stochastic unitary operator W+i on Hn+i by 


(4.7) 


W+l — Vn+lfi © Hn+1,1 © • • • © Hn+1,4"-1 

The eigenvalues of W+i are 1 (with multiplicity 4”'), j = 

0,1,..., 4”'—1. The corresponding unit eigenvectors are j = 0,1,..., d”— 

1, k = 0,1, 2, 3 given by fl4.7p . The operator W+i provides an intrinsic sym¬ 
metry on Hn+i generated by the coupling constants. We call W+i the cou¬ 
pling constant symmetry operator. In this framework, we can also dehne 
decoherence functionals and quantum measures in a standard way mn- 
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5 Dirac Operators 

In previous sections we have represented the c-causets in Vn by Xnj, n = 
0,1,..., i = 0,1,..., 4" — 1. We then denoted the standard orthonormal 
bases for the Hilbert spaces Hn by Xnj- In this section, we follow Dirac to 
dehne an energy-momentum operator on Hn- To accomplish this, it is useful 
to write j in its quartic representation 

3 = Jn-ld" ^ -|- jn- 24 "' ^ + ■ • • jl4 -|- Jo, ji ^ {0, 1, 2, 3} 

which as usual, we abbreviate 

j = 3n-ljn-2 ■ ■ • jiJo, ji e {0,1, 2, 3} (5.1) 

This representation describes the directions that a path turns when mov¬ 
ing from xofl to Xnj- For example, fl5.ip represents the path that turns in 
direction jn-i at a:o,o, then turns in direction jn- 2 , ■ ■ and hnally turns in 
direction jo before arriving at Xnj- If j has the form fl5.1l) and k G { 1 , 2 , 3}, 
dehne 

3 ~ 3n-l3n-2' ''3i3o 

where j' = ji if ji = k and j' = 0 if j) 7 ^ k. We also dehne j° = j. Notice 
that j° = j^ -I- j^ -I- j^. 

On Hn dehne the operators Qn^nj = and extend by linearity. 

We then have that 

iQir = iQir+{Qir+iQir 

Thinking of the Qn as coordinate operators, we dehne the energy-momentum 
operators Pn = VnQnyn-, A: = 0,1, 2, 3 on Hn- In dehne the indehnite inner 
product 

q-p = qopo - qiPi - q2P2 - qsPa 
and let (Ja;, A; = 0,1, 2, 3, be the Pauli matrices 



'1 

O' 


'0 

1' 


'0 -1 


[1 0 1 

do = 

0 

1 

, cn = 

1 

0 

, (T2 = 

i 0 

, 0-3 = 

0 

1 


Also dehne a = (cto, cti, (T2, <73) and = (P°, P^, P^, P^). We then have 

■ £„ = 

^ ' PS-PS -PS + iPS 
-PS-iPS PS + PS 
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We consider q_ ■ P„ as a self-adjoint operator on Hn 0 and call a ■ 
the discrete Weyl operator. To find the eigenvalues and eigenvectors of a ■ 
P„ notice that the eigenvectors of P^ are ynj = VnXnj with corresponding 
eigenvalues ^/j^. 

Theorem 5.1. Define finj = filj) and J = 0,1, • • •, 4*^ 

1, in Hn 0 by 



Then (j)n,j are eigenvectors of q_ ■ P„ with eigenvalue 0 and ifnj are eigenvec¬ 
tors of a ■ P_n with eigenvalues 2a/j°. 

Proof. A direct verification gives 



and 


- {P^ + + (P° + 

= (v0+ V?) (--P„‘ - iP^)y.., + (v? + »vT) (-P° + 



= 0 


As before, we have 

(p° - p^x^. + (-p,,^ + ^p^x^. 

= - 'iVf) - Vf) yn,j + [-Vf - 

= 2 (VTV - *VjV) VnJ = ‘2^/f filj 


Vn 
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and 


- {p^ + + (p° + 

= (-Vj^ - * V^) yn,j + (-Vj° - V7^) (V/ + Vf) yn,j 

= (-j^ - f - f - f - 2 VTv) yn,j 
= (-2j° - 2 VjV) yn,i = 2^/f 

We conclude that g_ ■ P_n4>n,j = 0 and q_ ■ P_n'4’n,j = 2^^ tjjnj, j = 
0,1,...,4^-1. ’ ’ ’ □ 

The discrete Dirac operator is the 4-diniensional extension of the discrete 
Weyl operator. We hrst dehne the 4x4 gamma matrices by 

<70 0 1 .To —CTfcl , TOO 

■'o=|_0 -ffoj’ = oj- *^ = 1.2.3 

Notice that jj are self-adjoint matrices, j = 0,1,2,3. The discrete Dirac 
operator is 

1 ■ Pn = loPn - llPl - 72^n “ bS^n 


po 

n 

0 

^p^ 

+ P?i 

0 

pO 

^Pl - p^ 

-^Pl 


-‘Pi - Pi 

_p0 

n 

0 

-^Pl + P?i 

‘Pi 

0 

_p0 

n 


We interpret the self-adjoint operator 7 ■ P„ as the empty-momentum oper- 
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ator. To present the eigenvalues and eigenvectors of 7 ■ we define 



{V2-l)^ 


0 


0 


(72-1) 

Ml = 

iVf 

, U2 = 

+ 7^ 


i\/¥ - 


-*77 




-iVf 


-iVf 



U 3 = 

0 

, M4 = 

(72-1)7? 


(72-1) 


0 


Theorem 5.2. The eigenvalues of 7 ■ P„ are — 1 ^ 2 ^ and j = 

0,1,..., 4” — 1. These eigenvalues have multiplicity two and their corre¬ 
sponding eigenvectors are upijng, K = 1,2, 3,4. 

Proof. By direct verification we have 


7 • Pn'^lVnJ 


{V2 - 1) f - f + - V^) 

(\/2 - 1) ^/f (^-i\/f^ - i\/ff 
(\/2- 1) ^/j0 ./f j - y/jO (^iy/j^- 


(V2-2)f 

0 

-i\/2ff 


UnJ 
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= U2yn,j 


7 • PnUsVnJ 



- (\/2 - l ) a / jV 




UnJ 


{iV¥+ V^) (-*V^+ ^/f) +f - {V2-l)f 



yn,j = V^U^ynJ 


j^+f+f - {V2-l)f 
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f +J^+f - (\/2 - 1 ) 


0 


Another way to express the results of Theorem 15.21 is the following. We 
can write 7o = <^0 <8) cr^ and 7fc = ® 0 ' 2 , k = 1.2.3. We also have that 

En ■ 1= Pn ® 70 - Pn ^ 'll - Pn ® 72 “ 73 

On Hn (8) For / G Hn, n G with v = {vq, Ui, V 2 , V 3 ) we have that 

3 

En ■ lf{Xn,j) ® -y = P^f{Xn,jhoV - ^ f {XnjhkV 

k=l 

= (V7° 7o - VE 7i - Vf 72 - Vf 73 ) V 



To hnd the eigenvalues and eigenvectors of P„ • 7, let / = ynj- Then 

E ■ XVnJ <^V = yn,j 
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where j = \/Pi \/Pj- Hence, we require solutions to j ■ 'yv = 

Xv. We conclude that the eigenvalues are —^/2p and with correspond¬ 
ing eigenvectors yn,j ® Uk, k = 1,2, 3,4. 


6 Multiverse Probabilities 


This section proposes specific values for the coupling constants. These values 
result in multiverse probabilities that predict a preponderance of pulsating 
universes. Although these universes are not exactly like the ideal pulsating 
universes considered in Section 3, they have similar properties. Whether 
these coupling constants provide an approximation to general relativity will 
be left for later studies. 

In the most elementary case, the coupling constants would be in¬ 
dependent of n,j so there would only be four constants c^, k = 0,1,2,3. 
Probably the simplest nontrivial values for these constants are m: 

c° = cos^ = -i cos 9 sin = - sin^ 


We can replace by 1 because this is an overall phase factor that does not 
affect probabilities. We can then write 

c^ j = cos^ 9, j j = -i cos 9 sin 9, = - sin^ 9 (6.1) 

where 9 G (—vr, vr). 

For j = 0,1,..., 4"' —1 we write j in its quartic representation fIS.ip . Then 
Xn,j is reached by the path {jn-i,jn- 2 , ■ ■ ■ ,ji,jo)', that is, the path 


(^ 0 , 0 ) • • • ) l+in-24" 2-|- 

For example 50 = 3 • 4^ -|- 0 • 4 -|- 2 so we have the quartic representation 
50 = 302. Hence, xs^so is reached by that path (3, 0, 2) which we can write 
as 

(^0,0) 3:1^3, X2,i2; ^3,50) 

The amplitude of Xnj becomes 


^ — Jn-lJn-2 Jn-3 


■ ■ ■r’° 

-2 n-lj„_l4" '^+jn-24X- 3-1- 
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For example, a(x 3 ^ 5 o) = Co gC° 3 C 2 2 - Now let j{k) be the number of ks in the 
quartic representation of j where k = {0,1, 2, 3}. If have the form fl6.ll) 
we have that 




For example, in this case we conclude that 


a(a^3,5o) = c^njcljclj = *cos^^sin^0 = |(sin20)3 

In a very simple model, let us assume that ^ have the form fl6.1l) and 
that 6 is moderately small, say 6* = 1/10 radians. Then 


cV = cos=(l/ 10 ) « (1 - jJj)" = ( 99 / 100)2 
<ij = c.lj = -!cos(l/ 10 )sin(l/ 10 ) « -j (l - j^) 
clj = -sin2(l/10) -1/100 
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-i/lO 


In the usual quantum formalism, the probability that the system is at x. 


n,J 


IS 




g |2i(0) I 1 l2j(l) I 2 |2i(2) I 3 121(3) 


n,J I 


n,j I 


n,j I 


n,j I 


Suppose n is large. Then among the c-causets Xnj, we have that x„^g has 
the largest probability and has essentially zero probability. However, 

there are a large number of other c-causets (paths). In particular, there are 
many c-causets for which j(0) is large, j(3) = 0 and j(l) and/or j(2) are small 
but nonzero. Taken together, these c-causets have a dominate probability. 
Such c-causets have approximately the geometry of a pulsating universe as 
considered in Section 3. For example, suppose n = 20 and j(0) = 17, j(l) = 
1, i(2) = 2, j(3) = 0 with 


j = 20000200000001000000 


so that j = 2 ■ 4^®-1-2 • Then Xnj has the shell sequence (1,10,17,13). 

To obtain a more realistic model, we propose that 6nj is a function of n 
that provides a more periodic behavior. This would result in a closer approx¬ 
imation to a pulsating universe as considered in Section 3. In particular, we 
suggest that cj; ^ have the form (16.Ih with 


e 


n 


1 

ln(n -|- 1) 


cos [ln(n -|- 1)] 
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We would then have 


Cn,i = cos 2 ( 0 „) = cos^ cos [ln(n + 1 )] 

Cn,i = sin cos [ln(u + 1)] 

— 1 

Upon graphing these functions one will see a kind of periodic behavior with 
increasing periods. 
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